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Abstract 

We present explicit reconstruction algorithms for fully anisotropic unknown elasticity 
tensors from knowledge of a finite number of internal displacement fields, with applications 
to transient elastography. Under certain rank-maximality assumptions satified by the strain 
fields, explicit algebraic reconstruction formulas are provided. A discussion ensues on how to 
fulfill these assumptions, describing the range of validity of the approach. We also show how 
the general method can be applied to more specific cases such as the transversely isotropic 
one. 


1 Introduction 

We consider the reconstruction of a fully anisotropic elasticity tensor C = from 

knowledge of a finite number of displacement fields solutions of the system of linear 

elasticity 


V • (G : (Vu + (Vu)^)) = 0 (A), = g (prescribed). (1) 

Applications for such a theory include the medical imaging modality called elastography. 
Elastography is concerned with the reconstruction of the elastic properties in biological tissues. It 
has been observed experimentally that certain biological tissues (e.g., muscle fiber [19], or white 
matter inside the brain [18]) display anisotropic mechanical properties, and the present article 
aims at giving access to these anisotropic features. The present approach to elastography consists 
of two steps. A first step is the reconstruction of the internal elastic displacements, which are 
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accurately modeled as solutions of the above system of linear elasticity. Two different methods 
are typically used to reconstruct such displacements. One such method, called Ultrasound 
Elastography, consists of probing the elastic displacements by sound waves. We refer the reader 
to [20, 32, 33, 34] for more details on the method. A recent analysis of such reconstructions from 
ultrasound measurements is performed in [11]. A second method, called Magnetic Resonance 
Elastrography, leverages the displacement of protons by propagating elastic signals to reconstruct 
the elastic waves by magnetic resonance imaging [21, 22, 30]. 

In this paper, we assume the first step done and consider the quantitative reconstruction of 
elastic coefficients from knowledge of such elastic displacements. For references on this quanti¬ 
tative step of elastography, we refer the reader to [15, 6, 27]. Here, we generalize earlier work 
performed for scalar second-order equations [13, 12] to the case of systems, and we generalize 
recent work on the Lame system [6, 24] to the fully anisotropic setting. Other hybrid inversions 
for elliptic systems can be found in the case of Maxell’s system in [8, 14] and in a more general 
framework in [5]. 

The approach consists in deriving explicit algebraic formulas reconstructing the unknown 
parameters locally, based on hypotheses of linear independence (or rank maximality) of func¬ 
tionals of the measurements and their partial derivatives. A discussion then follows on what 
regularity or property (e.g., the Runge approximation property) is required a priori on the un¬ 
known parameters so that the hypotheses of reconstructibility may be fulfilled. The idea of 
constructing local solutions fulfilling certain maximality conditions, which are then controlled 
from the boundary of the domain via Runge approximation, was also used in the context of 
reconstruction of conductivity tensors from knowledge of so-called power density functionals [29] 
or current density functionals [9]. 

The above problem of reconstruction of coefficients in partial differential equations from 
knowledge of internal functionals is sometimes referred to as a hybrid inverse problems; see 
[2, 3, 4, 5, 23, 35] for reference for hybrid inverse problems in other imaging modalities. Note that 
we do not consider here the inverse boundary elasticity problem, for which Lipschitz stability 
estimates may be obtained only when the Lame coefficients are piece-wise constant; see for 
instance [16] for a recent reference on such a topic not covered here. 

We now give the main results of the paper in the next section and give an outline of the 
remainder of the paper there. 

2 Main results 

Preliminary notation and definitions. In what follows, we denote by M 3 (M) the vector 
space of 3 X 3 real matrices with inner product A : B := ti (AB'^) = with 

respect to which we recall the orthogonal decomposition M 3 (M) = 53 (M) © A 3 (M), where the 
first (second) summand denotes (skew-)symmetric matrices. 

Let us fix A C a bounded domain with smooth boundary for the remainder of the paper. 
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An elasticity tensor C is a fourth-order tensor satisfying the following symmetries 

Cijkl Cjikh Cijkl Cijiki Cijkl Ckliji 1 ^ ^ 3, 

characterized by 21 independent components (instead of 81), where the latter symmetry corre¬ 
sponds to the assumption that C is hyperelastic. We assume below that C is uniformly pointwise 
stable over X [26, Ch. 6, Def. 1.9] in the sense that there is a k > 0 such that 

: C'(x) : e > K € : e, VxGX, VeG 5'3(M). (2) 

For C = {Cijki}ij,k,i an elasticity tensor and e € S'slM), we denote C : e := {Ylk,i Cijkieki}i,j- 
In some sections below, the tensor C will be represented in the non-tensorial Voigt notation, 
an 56(M)-valued function c = {c„^}i<„^^<g, where the correspondence of elements Cq,^ in terms 
of the coefficients Cijki is obtained via the double index mapping lli—)-l, 22i—)-2, 33i—>-3, 
23,32 1 -^ 4, 13,31 i—)• 5 and 12,21 i-)- 6 (for instance, cn = Cim and so on). Hooke’s law, 
relating stress a to strain tensors e via the relation a = C : e, now reads in Voigt notation 
cry = c ey, where we define 

ey = (eii, £ 22 , £33, 2623, 2631,2ei2)^, = (<7ii, 0 - 22 , 0 - 33 , £^23, <731, £^ 12 )^- (3) 

Most often, we will drop the subscript “V” below as the context will tell us naturally what 
representation to pick. For (6^^\ ..., 6^®)) in 53(M), we define below 

det(e«),...,el«l) = det{£l;>,...,e®). (4) 

y R<3 

A crucial fact for further derivations is the following. 

Lemma 2.1. Under assumption (2), there exists k' > 0 such that 

detc(x) > k', V X € A. (5) 

Main results. Provided linear independence conditions that can be satisfied locally and 
checked directly on the available measurements, we first provide an explicit reconstruction al¬ 
gorithm for fully anisotropic elasticity tensors. Here and below, a typical displacement field 
is denoted u : A ^ with corresponding strain tensor 6 : A —>■ 53 (M) with components 
Cij = ^{diUj + djUi), or, in coordinate-free notation, e = ^(Vu -|- (Vu)^). 

We now formulate our main hypotheses in order to set up our reconstruction procedure. These 
hypotheses, based on algebraic redundancies of various elasticity solutions, force the unknown 
tensor to lie on the orthogonal of a space generated by a rich enough set of data. The hypotheses 
below formulate how some functionals of the data set can be made to generate a hyperplane of 
S'6(M), a normal of which can be explicitely constructed and proved to be proportional to C, 
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the proportionality factor being reconstructed at the end. As seen below, fulfilling this agenda 
requires 6 + solutions, where N > and d denotes the number of scalar components of 
C (the factor 3 accounts for the fact that each additional solution after the sixth one provides 
three orthogonality constraints on C). 

Hypothesis 2.2. Let Li C X. 

A. There exist 6 solutions of (1) whose strain tensors form a basis of S^fR) at 

every x € fl. This condition can he summarized as (dety is defined in (A)) 

inf det(e^^^(x),..., e^®^(x)) > cq > 0, for some constant cq. (6) 

xGO V 

B. Assuming A fulfilled, there exists N additional solutions giving rise to a 

family M of 3N matrices whose expressions are explicit in terms of {e^^\daC^^\ 1 < a < 
3, 1 < j < 6 +A^} (see (12) for detail), and such that they span a hyperplane of SefR) at 
every x € 12. This condition can he summarized as 

inf ^ M{M') : > ci > 0, for some constant ci, (7) 

M'CM, #M'=20 

where M is an operator generalizing the cross-product, defined in (13). 

Remark 2.3. It should be noted that these hypotheses are stable under smooth perturbation 
of the boundary conditions generating the displacement fields , ■ ■ ■, . This is because 

Hypotheses 2.2. A-B are expressed in terms of functionals which depend polynomially on the 
components of displacement fields and their derivatives up to second order, which are in turn 
continuous functionals of their boundary conditions (see, e.g.. Sec. 5.2.1 for appropriate topolo¬ 
gies). 

Remark 2.4. In Hypothesis 2.2.B, the number N depends on the type of isotropy of the tensor 

C. In the most general, 21-parameter case, spanning a hyperplane of Sq{R) with 3N additional 
constraints suggests N >7, and the total number of displacement fields needed is then 6 + 7 = 13. 

On an open subset 12 C A where Hypotheses 2.2. A-B hold, we then derive an explicit 
reconstruction algorithm, reconstructing C over 12. This is done as follows: in Voigt notation, 
decompose the S'6(M)-valued function C as the product of a scalar function r times a normalized 
anisotropic structure C such that C = tC {C is such that its Voigt counterpart c has determi¬ 
nant 1, so r = (detc)®). Algebraic manipulations allow us to obtain pointwise orthogonality 
constraints which, under hypotheses (6)-(7), are numerous enough to reconstruct C pointwise. 
Once C is reconstructed, an equation for Vlogr is derived over 12, leading to the reconstruction 
of r by solving either a transport or a Poisson equation. Finally, once C = tC is known, ad¬ 
ditional stability is recovered on certain components of C by deriving equations reconstructing 
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the third-order tensor div C (whose components may be written as diCijki, 1 < j,k, I < 3). 
Such components, although linear combination of first-order derivatives of the components of C, 
are reconstructed with the same stability as C itself. Such stability improvements have already 
been observed in e.g. [10, 9]. 

The approach just described yields unique and stable reconstructions over in the sense of 
the following theorem. 

Theorem 2.5. Suppose that over some open set Q G X, hypotheses 2.2.A-B hold for two 
families of displacement fields and corresponding to elasticity tensors C 

and C . Then C and C can each he uniquely reconstructed over Q from knowledge of their 
corresponding solutions, with the following stability estimate for every integer p > 0 

N+6 

\\C — C'^||vKP"°°(o) + II C — div C''||wp.°o(o) < K ^ — e (8) 

j=i 

Remark 2.6. Note that if the elasticity tensor is split into the product of an unknown scalar 
function r times a known anisotropic tensor, the stability of the problem of reconstructing r 
from fields is better-posed, i.e. involves the loss of one derivative instead of two, according 
to the statement 

N+6 

Ik “ k||wp+i.°°(r2) < K ^2 ^ ^'^^llvKp+i'°°(0)- (9) 

1=1 

In the context of this reconstruction approach, an elasticity tensor is then reconstructible 
on some given open set Q (or globally on X) if there exist displacement field solutions of (1) 
fulfilling Hypotheses 2.2.A-B throughout fl. In this context, classifying reconstructible elasticity 
tensors then consists in finding out in what situations one can construct such displacement fields. 
Two such situations are described in the following theorem: (i) a fully anisotropic tensor that 
is close to constant is globally reconstructible from well-chosen displacement fields whose traces 
are explicitly given, and (ii) an elasticity tensor with smooth enough components which satisfies 
the Runge approximation property (see Section 5.2.3, in particular, Eq. (39)) is, in principle, 
locally reconstructible from knowledge of its displacement fields. 

Theorem 2.7 (Reconstructibility of elasticity tensors). In either of the following cases, there 
exists a non-empty open set of smooth enough boundary conditions generating displacement fields 
characterizing C uniquely and stably in the sense of Theorem 2.5. 

(i) C is -close to constant. 

(ii) C is smooth (at least of class C^) and satisfies the Runge approximation property. 
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An application of particular interest is the reconstruction of transversely isotropic (TI) elas¬ 
ticity tensors. Such tensors have, as of yet, the highest type of anisotropy for which the Runge 
approximation has been proved (see [31]), so that they are covered in case (ii) of Theorem 2.7. 
We also explain in Section 5.1 how to generate explicit boundary conditions reconstructing a 
constant TI tensor which, by case (i) of Theorem 2.7, can be used to reconstruct a near-constant 
TI tensor. 

Outline. The rest of the paper is organized as follows. Section 3 covers the proof of Lemma 

2.1. Section 4 covers the reconstruction procedure as well as its stability property (proof of 
Theorem 2.5). Section 5 covers the two ways described in Theorem 2.7 to fulfill Hypotheses 

2.2. A-B for certain classes of tensors, thereby establishing their unique and stable (in the sense 
of Theorem 2.5) reconstructibility from displacement fields. 


3 Proof of Lemma 2.1 


Suppose the elasticity tensor C satisfies (2). In the study of elastic eigentensors in [28], it 
is established that C can have at most six distinct eigencouples (e, A) such that the relation 
C : e = Ae is satisfied. The symmetries of C and hypothesis (2) imply that the A’s are real and 
that all of them satisfy A > k. It is also mentioned in [28] that the eigenvalues A are precisely 
the eigenvalues of the 56(M)-valued representation of C in the form 


c' = {2^' Cjj}i<ij<6, where x{i) 


0 if i = 1,2,3, 
1 if i = 4, 5,6, 


and where Cij denotes the Voigt representation of C presented in the introduction, 
immediate that detc' = 8detc. Moreover, as mentioned above, all eigenvalues of c' 
eigenvalues of C and therefore satisfy the estimate A > k. Therefore 

, 1 , 
det c = - det c > —, 


It is then 
match the 


hence (5) holds with k' = 


4 Reconstruction algorithm and its stability 

As seen in Theorem 2.5, the left hand side of (8) contains two terms whose product forms the 
elasticity tensor C, reconstructed with stability estimates in different norms: we decompose C 
into the product tC, where C contains the rescaled anisotropic structure of C, defined by a 
normalizing condition (specifically, if c describes C in Voigt notation, then det c = 1) and r is 
the remaining scalar factor. That this is possible comes from estimate (5), which states that 
detc is uniformly bounded away from zero throughout X. The next two sections focus on the 
successive reconstruction of C, then r upon assuming that Hypotheses 2.2.A-B holds. 
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4.1 Reconstruction of the anisotropy 

In this section, we will use the Voigt notation, so that strain tensors will be represented as Re¬ 
valued funtions e = ey. In this notation, the elasticity tensor becomes an S'6(R)-valued function, 
characterized by 21 scalar functions c = {cap}i<a,p<Q satisfying Cq,^ = c^q, and where Hooke’s 
law is expressed as a regular matrix-vector product ay = c ey, with ((Ty,ey) as in (3). This 
makes the problem tractable via algebraic manipulations on matrices instead of 4-tensors. 
Using the Voigt notation, the elasticity system (1) takes the form 

/ di 0 0 0 <93 <92 \ 

Dy • (c e) = 0, Dy ■.= { 0 (92 0 03 0 (9i j , (10) 

\ 0 0 (93 02 <91 0 / 

and where, for a Me(R)-valued function A and a scalar function /, we have the identity 

Dy{fA) = {Dyf)-A + fDyA. (11) 

As the reconstruction approach is local (even pointwise for the anisotropic part), we assume 
to have 6 elasticity solutions whose corresponding strain tensors {ey^(x)}i<j<6 form 

a basis of R® for every x of some subdomain H C A. Any additional solution with p > 7 is 

such that for x € 11, decomposes uniquely into the basis above as 


= where ■; .f (x)) ^ 

J = 1 


dety(e(^) (x),..., e(®) (x)) 


Plugging this equality into the elasticity equation, we obtain that 

6 6 

0 = By (ce^P^) = '^Dy {npjCe^^\xo)) = '^{Dy^pj) ■ -7 [ipjDy ■ (ce^^^), 

i=i i=i 

which, since Dy ■ (ce^-^^) = 0, implies that 

6 

'^{DyUpj) ■ = 0 . 

1=1 


This last equation can be seen as three scalar orthogonality constraints on the tensor c in the 
inner product structure of S'6(R) which we denote by A : H := tr [AB'^) = j=i where 

the matrices that are othogonal to C are directly known from available measurements. The last 
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equation is equivalent to 


c : = c : = c : = 0, where 

6 

jyjlp),! ^ ^(9i/ipj 0 0 0 ds^ipj d2fipj) (8) 

i=i 

6 

^,^(p),2 ^ Q Q (g) gO)^ 

J=1 

6 

m(^)’3 = ^(0 0 ds^ipj d2iJ.pj difipj 0) (g) 
i=i 


( 12 ) 


Note that since c is orthogonal to ^6(1^)) one could replace the matrices with their 

symmetrized versions. Notice that the components of these matrices are first partial derivatives 
of rational functions of strain tensors. 

If enough linear constraints of the form (12) are available from a rich enough set of measure¬ 
ments, that is to say, if enough such matrices are available and form a hyperplane of 56(M) at xq, 
then the tensor c(xo), constrained to be perpendicular to this hyperplane, will be determined 
up to a multiplicative constant. The reconstruction can be done via a generalization of the 
cross-product, as used for instance in [29] for similar purposes in the context of the conductivity 
equation. Define {mj }fLi a basis of Sq{M) and given M = {Mj}^2=i C 5'6(M), define the oper¬ 
ator J\f : S'6(M)^^ —>■ 5'6(M) by expanding the formal determinant below with respect to its last 
row 


Ml : mi 


N{M) := 


1 

det(mi, • • • ,m 2 i) 


M20 : mi 
mi 


Ml : m2i 


M20 : m 2 i 
m2i 


(13) 


AA is a 20-linear, alternating map that does not depend on the choice of basis for 56(M). Af{M) 
is a vector that is normal to the hyperplane spanned by M when M is linearly independent, zero 
otherwise. In particular, if M is a family of matrices known to be orthogonal to a given matrix 
m', then M{M) is either zero (if dim span M < 20) or proportional to m' (if dim span M = 20). 


In light of this last comment, and assuming that a rich enough set of solutions of the elasticity 
system (1) gives rise to a family of matrices M of the form (12), with cardinality greater than 
20 and spanning a hyperplane of 56(M) at a given point xq G X, then for any given 20-uple 
M' C M, M{M') is either zero or proportional to c(xo). Normalizing c(xo) = det(c(xo))~^c(xo) 
and enforcing the condition cn > 0 (this is because in the proof of Lemma 2.1, we notice that 





the top-left 3x3 block of c matches that of c' which is a symmetric positive definite matrix, so 
one must have cn > 0, C22 > 0 and C33 > 0), then we have the equality 

{±)m'N{M') = (det(Ar(M')))^c(xo), 

for every 20-uple M' C M, where (±)m' is the sign of the top-left entry of This equation 

is either trivial when M' is linearly dependent, or reconstructs c(xo) otherwise. Condition 7 
ensures that at least one subfamily M' is linearly independent, so we sum the last equation over 
all subfamilies and arrive at the formula 

c(x„)=( (det(V(Af')))4) (±)„W(Af'). ( 14 ) 

\M'CM, #M'=20 / M'CM, #M'=20 

Remark 4.1. The expliciteness of formula (14) is interesting in its own right and makes the 
stability of the problem straighforward to assess. On the other hand, the computation of 20 x 20 
determinants can be expensive, and methods constructing a normal to a hyperplane without 
implementing (14) might reveal more practical. For a 20-uple M, an example of a potentially 
faster method for finding a scalar multiple of M{M) is to form a 20 x 21 matrix whose rows are 
the elements of M, and to find a vector in the nullspace of that matrix via Gaussian elimination. 

Now that the anisotropic structure C (or equivalently, c) is reconstructed, we now explain 
how the multiplicative scalar r can be reconstructed via a standard transport equation. 

4.2 Reconstruction of the scalar factor r 

We now switch back to 4-tensor notation. Plugging the decomposition C = tC into the elasticity 
equation, where C is assumed to be known from the previous step, we obtain, for each elasticity 
solution considered, 


{C : e)VlogT = —div {C : e), (15) 

where (C : e) G 5 ' 3 (ffi) is not necessarily invertible. However, under the assumption that 

is a basis of S' 3 (M), by virtue of (5), so is {C : Therefore, the identity tensor I 3 

decomposes into this basis by means of some functions //j(x) such that 

6 _ 

l3 = j;/r,(x)C':e(^). 

1=1 

If we denote Dij := (C : e^*^) : (C : for 1 < < 6 , the S' 6 (M)-valued function D = {Dij} 

is known and invertible (as the Grammian matrix of the basis {C : and the entries of 
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its inverse are denoted by In this case, the functions take the explicit form 

6 _ ® _ 

^, (x) = I 3 : (C : e«) = tr (C : 6 »). 

i=l i=l 

Now taking a linear combination of (15) weighted by the functions /Xj, we obtain 

6 _ 

V log T = — E ^j(x) div {C : (16) 

i=i 

The right hand side of this equation is completely known and r can therefore be reconstructed 
via either {i) integration of (16), or, after taking divergence, {ii) solving a Poisson equation on 
the domain Q with known Neumann boundary conditions. As discussed in [7], the approach (ii) 
may be more robust to noise than the former, as the integration of ordinary differential equations 
in the presence of noisy measurements may strongly depend on the choice of integration path. 
Moreover, taking divergence of (16) has the advantage of naturally removing the curl part of 
noisy right-hand sides in (16). 

4.3 Reconstruction of the tensor div C 

Now that both C and r are reconstructed, we explain how to gain stability on the reconstruction 
of the third-order tensor div C = diCijki ej ( 8 ) (g) e;. Note that this tensor is symmetric in the 

pair of indices (k, 1). Now the elasticity equation can be rewritten as 

{diCijki^^ki T 0, 1 E J E 3, 


or in contracted notation. 


(div C)j.. : e = -Cijkidieki, 1 < j < 3. (17) 

For every 1 < j < 3, (div C)j.. can be seen as an 53(M)-valued function, and as such, using again 
the assumption that is a basis of 5'3(M), upon defining Eij = e*-*) : for 1 < i,j < 6, 

the S'6(M)-valued function E = {E^j} is uniformly invertible over and we denote by E'^^ the 
components of its inverse. With such notation, Cramer’s rule in S' 3 (M) then reads 

(div C)j.. = EP‘1 {{div C)j.. : 1 < J < 3, 

and combining this with (17), we obtain a reconstruction formula for the tensor div C 

(div C)y. = -EP^{Cijki dte^^) 1 < J < 3. (18) 
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4.4 Proof of Theorem 2.5 

Uniqueness. Equations (14), (16) and (18) give us explicit algorithms to reconstruct c, r and 
div C explicitely, with the only indeterminacy that r is defined np to a multiplicative constant. 
This is becanse multiplying the elasticity system by a constant changes neither the eqnation, nor 
its solutions. We can remove this indeterminacy by assuming that r is known at a (boundary) 
point or by fixing its value at a point. The constructive nature of the approach gives uniqueness. 

Stability on C. The functional reconstructing c is a rational function of strain tensors and 
their first partial derivatives. If two sets of displacement fields and {u both 

satisfy hypotheses 2.2.A-B over some Q <Z X, the denominator of the right hand side of (14) 
never vanishes, and its rational expression in terms of measurement components yields the local 
estimate 

_ _ e+N 

l|C - C'lUoo(f,) <kY, (19) 

i=i 

for some constant K > 0. More generally, if the strain tensors are in IT^'’'^’°°(n) for some integer 
p > 0, the reconstruction formula (14) can be differentiated p times to yield stability estimates 
in smoother norms of the form 

_ _ 6+Ar 

\\C - C'\\wp,°°{n) < K P > 0, (20) 

where the constant K grows polynomially in terms of the maximum lT^^^’°°(fI) norm of the 

strain tensors and the inverse of min(co, ci, Cq, c'^) > 0, where (cq, ci, Cq, c'j^) are the constants 
defined in (6)-(7) corresponding to each system of measnrements and {u 

Stability on the scalar factor r. Equation (16) takes the form Vlogr = F{C,u) with F 
a rational function of the components of C and its first derivatives, and of strain tensors and 
their first partial derivatives, one deduces 

6+Af ^ ^ 

Ik ~ ^ f lliyp+i’°°(o) + K"\\C — C'kwp+i’°°(o)) P > 0. (21) 

i=i 

We see that when C is known a priori, only one derivative is lost from the measurements to 
the quantity r. However, when considering joint reconstructions of (C,r), we see that errors 
\\C — C'l|i 4 /p+i.oo(Q) are controlled by measurement errors in — u norm, so 

that the second term in the right hand side of (21) implies a loss of two derivatives from the 
measurements to the quantity r. This explains Remark 2.6. 
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Stability on div C. The right-hand-side of (18) is a linear functional in the components of C, 
rational in the components of the measurements and its partial derivatives up to second order 
so that, as before, the stability on div C is of the form 

_ 6-l-Af 

||divC-divC||ivp.oc <K'Y, P > 0- (22) 

i=i 

This ends the discussion on stability and the proof of Theorem 2.5. 

5 Reconstructible tensors - Fulfilling Hypotheses 2.2.A-B 

5.1 The constant coefficient problem 

We now show that, in the same way that harmonic polynomials up to second order uniquely 
characterize a constant diffusion tensor in the scalar case, a constant elasticity tensor C is 
uniquely characterized by elasticity solutions whose components are polynomials up to second 
order. 

Polynomial displacement fields. Assume that C is constant and let us denote by c* the 
i-th row (or column) of c (i.e., C in Voigt notation). 

ft is straightforward to see that any displacement held u with linear components is a solution 
to (1). These solutions in fact allow us to construct a basis of strain tensors, i.e. such that 
..., is the natural basis of for every x € A. This can be achieved by considering the 
following solutions 

uW(x) = (x,0,0)^, u(2)(x) = (0,y,0)^, u(3)(x) = (0,0, z)^, 

uW(x) = i(0,z,y)^, u(5)(x) = i(z,0,x)'^, u(6)(x) = i(y,x,0)^. 

Second, a displacement held with quadratic components, of the form 

u(x) = Qx • Px, ^x • Qx, ix • Px^ , P,Q,R£ 53(M), 

has a strain tensor, in Voigt notation, taking the form e(x) = xVi + yV 2 -|- zV^, where we have 
dehned 

Vi ■= (Pn, Q21, R31, Q31 + R21, P13 + R11, P12 + Q11)"’", 

V2 := (P12, Q22, R32, Q32 + R22, P23 + R21, P22 + Q2i)^, 

V3 := (Pi3, Q23, R33, Q33 + R23, P33+R31, P 32 + Qi 3 )^- 


(23) 


(24) 


12 



The correspondence {P,Q,R) i->- (Vi,V 2 ,V 3 ) is bijective, with inverse {Vij denotes the j-th 
component of V^) 



1 

i-H 

1 _ 


■ Vhe - ^21 ^12 UVse + Via - V25) ' 

p = 

V26 — V12 2(^25 + ^36 — bi 4 ) 

, Q = 

V22 

^32 


. sym • V35 - Vi3 


sym 

1 

CO 

1 _ 


R = 


^15 - ^31 i(^14 + ^25 - ^36) ^13 ' 

V24 — V32 V23 , 

sym • 1/33 _ 


(25) 


so that prescribing one or the other is equivalent. Then the stress tensor is expressed as 


a = c e = X c Vi + y c V 2 + z c V 5 . 

With this expression, u is an elasticity solution if and only if the three scalar conditions below 
are fulfilled 


Cl • Vi + C6 • 1/2 + C5 • VJj = 0, 

C6 • Vh + C2 • V2 + C4 • V3 = 0, ( 26 ) 

C5 • Vi + C4 • V2 + C3 • V3 = 0 . 

This essentially leaves us with a 18-3 = 15-parameter family of quadratic solutions to the 
constant-coefficient elasticity system, and we now show that these solutions suffice to char¬ 
acterize C uniquely. Conditions (26) can be written in the form (12) where, defining the 

natural basis of M®, the matrices M^, M^, take the form 

= Vi (g) ei -h 1/^2 <8) ce -h V3 (g) 65, 

= Ti (g) eg -h 1^2 (g) 62 -h V3 (g) 64, ( 27 ) 

= Vi (g) 65 -h 1/2 <8> e4 -h V3 (g) 63. 

We will show in the next paragraph that, for various choices of Vi,V 2 ,V 3 satisfying (26), the 
corresponding matrices (27) span the hyperplane {c}-*- in M6(M), thereby imposing enough 
orthogonality conditions on c to determine it uniquely up to a multiplicative constant. 

Rank maximality of quadratic displacement fields. Since c is invertible, let us denote c* 
the i-th row (or column) of c“^, so that Cj • ct = Sij. Since we can write c as c = Y^j=i ej, 
using the identity 

{U 0 V) : {S (g)T) = {U-S){V -T), U,V,S,T€R^, 
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we can show that the space {c}-*-, regarded as a hyperplane of M 6 (M) of dimension 35, is spanned 
by the following family: 

{c}-*-= span {c* (g) Cj, 1 < i,j < 6 , i ^ j, c* (g e, — c*^;^ (g Cj+i, 1 < i < 5} . (28) 

Next, proceeding by exhaustion, we construct quadratic solutions giving rise to matrices 
Af^, spanning the family (28). 

• We first consider quadratic displacement helds satisfying (26), and such that two vectors 
among {Vi, V 2 , V 3 } vanish identically. 

- If Vi = V 2 = 0 , conditions (26) read Vs-Ljcs, C 4 , C 5 }, so that V 3 € span {c^, c^, Cg}, 
thus leading to matrices in the set {c* g) ej, i = 3,4, 5, j = 1, 2, 6 }. 

— If Vi = V 3 = 0 , conditions (26) read V 2 -L{c 2 , C 4 , cg}, so that V 2 € span {c^, Cg, Cg}, 

thus leading to matrices in the set {c* g) ej, i = 1,3, 5, j = 2,4, 6 }. 

— If V 2 = V 3 = 0 , conditions (26) read Vi_L{ci, C 5 , cg}, so that Vi € span {c^, Cg, C 4 }, 

thus leading to matrices in the set {c* g) e^, i = 2,3,4, j = 1, 5, 6 }. 

• Secondly, we consider quadratic displacement fields satisfying (26) and such that one vector 
among {Vi, V 2 , V 3 } vanishes identically. 

— If Vi = 0 , conditions (26) read cg • V2 + cg • V3 = C2 • 1^2 + C4 • V3 = C4 • V2 + C3 • V3 = 0 . 
One way to achieve this is by writing, for some free parameters a, (5, 7 , 

1^2 = ctCg + / 3 c 2 + 7 C 4 , V 3 = —acg — fjc\ — 7 Cg. 

The matrices M^, Af^, M'^ thus constructed take the form 

= a(cg g eg - C 5 g) eg) + / 3 (c 2 g eg - C 4 g eg) + 7(04 g eg - Cg g eg), 

Af^ = Q;(cg g e2 - C5 g e4) + / 3 (c 2 g e2 - C4 g e4) + 7(04 g e2 - Cg g e4), 

A/^ = a(cg g e4 - C5 g es) + / 3 (c 2 g e4 - C4 g es) + 7(04 g e4 - Cg g e^). 

— If V2 = 0, conditions (26) read ci • Vi + cg • V3 = cg • Ih + C4 • V3 = cg • Vi + C3 • V3 = 0. 
One way to achieve this is by writing, for some free parameters a, (5, 7 , 

Vi = ac^ + /Scg + 7C5, 1/3 = -acg - /3c4 - 7C3. 

The matrices M^, M^, thus constructed take the form 

= a{c*[ g ei - Cg g eg) + /^(cg g ei - C 4 g eg) + 7(05 g ei - Cg g eg), 

Af^ = q;(c 4 g eg - Cg g e 4 ) + /3(cg g eg - C 4 g e 4 ) + 7 (c 5 g eg - Cg g e 4 ), 

= a(c 4 g eg - Cg g e 3 ) + /3(cg g eg - C 4 g e 3 ) + 7 (c 5 g eg - Cg g e 3 ). 


14 


Examining the five cases considered, we see that the set (28) can be spanned by matrices 
generated by either of the above cases. This concludes the discussion. 

Remark 5.1. The number of solutions required here is not sharp. In fact, it would he enough 
to construct a hyperplane (of dimension 20 instead of 35) of Sq{E.) using well-chosen solutions, 
though it is not necessarily straightforward to find out which subfamily of (28) of cardinality 20 
spans the orthogonal space to {c} in 56(M). 

The transversely isotropic case. We treat here a particular example of anisotropy, where 
the number of unknowns is reduced to 5, and we show how the general method can be adapted. 
Assuming that the 03 (or “z”) direction is the constant direction of isotropy, we decompose a 


transversely isotropic tensor, in Voigt 

notation, 

as 




a 

b 

c 

0 

0 

0 


h 

a 

c 

0 

0 

0 


c 

c 

d 

0 

0 

0 

c = 

0 

0 

0 

e 

0 

0 


0 

0 

0 

0 

e 

0 


_ 0 

0 

0 

0 

0 

a—b 

2 


Locally, four well-chosen orthogonality constraints are enough to locate c(x) in the five¬ 
dimensional space describing it. As each displacement field gives rise to 3 orthogonality con¬ 
strains, we expect that two well-chosen displacement fields in addition to the six ones 

forming a basis of strain tensors should suffice. 

As in the fully anisotropic case, the family of linear displacement fields (23) forms a basis 
of strain tensors of elasticity solutions, and we now aim at finding two additional solutions in 
the form of well-chosen quadratic polynomial displacement fields. Since the unknowns are now 
the scalars (a, 6 , c, d, e), we can rewrite the conditions (26) for a quadratic displacement field to 
be an elasticity solution as orthogonality constraints on the vector (a, 6 , c, d, e). In this set of 
variables, a displacement field of the form (24) solves the system of elasticity if and only if 


a 


1^11 + ^^26 

Ih2 - ^V 2 e 

Ids 

0 

^35 


b 


■ 0 ■ 

kVw + ^22 

-kVl 6 + V 2 l 

1^23 

0 

1^34 


c 

= 

0 

0 

0 

V31 - 1 - V32 

1^33 

Vlb + V24 


d 


0 


e 


The two additional quadratic solutions can be chosen as follows: 
• is constructed so that 


■ Ifil + ^^26 

Ih2 - ^V2e 

Vl3 

0 

^35 


—b a 

0 

0 

0 ■ 

IViG + V22 

+ V 21 

V 23 

0 

V34 

= 

0 -c 

b 

0 

0 

0 

0 

V3I + ^32 

^33 

hl5 + ^24 


1 

0 

0 

-d 

c 

0 
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all other coefficients Vij being set to zero. One possible solution of this is 


(Vl 2 ,Vll,^ 23 ,V^ 33 ,^ 21 ,Vl 5 l) = {a,-b,b,C, -C, -d), 

(all other coefficients set to zero) which, upon using (25), yields the displacement field 
= {—bx^ — 2 cxy — ay^ — 2 dxz, cx^ + 2 axy — bz‘^, dx^ + 2 byz + cz^). 

• is constructed so that 


■ fil + ^V 26 

P12 - 1^26 

Pl 3 

0 

P35 


■ 0 

0 

0 

0 

0 ■ 

IE16 + P22 

-iEl 6 + P 21 

E23 

0 

E34 

= 

0 

0 

0 

0 

0 

0 

0 

V3I + V32 

P33 

bi 5 + P24 


0 

0 

0 

—e d 


all other coefficients Vij being set to zero. One possible solution is 

(^ 33 ,^ 15 ) = (-e,d), 

(all other coefficients being set to zero) which, upon using (25), yields the displacement 
field 


= ( 0 , 0 , - ez^). 

and were constructed to be solutions of the system of elasticity, and the orthogonality 
constraints they generate saturate a hyperplane of Satisfying such orthogonality constraints, 
(a, 6 , c, d, e) are uniquely determined up to a constant, which is in turn determined following the 
approach to reconstruct r in the general case (see Sec. 4.2). 

5.2 Proof of Theorem 2.7 

5.2.1 Preliminaries: forward theory and interior regularity 

Here and below, if u denotes a vector-valued function, we will loosely write u € H^(X) to mean 
that u has components in the Sobolev space H^{X). In order to prove Theorem 2.7, we will use 
the following theorem (see e.g. [26, Chap. 6 , Th. 1.11]): if C is smooth, uniformly pointwise 
stable (as defined in Eq. (2)) and dX is smooth, then for every / G L?‘{X) there exists a unique 
solution u G H‘^{X) to the problem 

V-(CVu) = f (X), u|ax = 0. (29) 

If f G H^{X) then u G H^~^‘^{X) for s > 0 with an estimate of the form 

l|u||_f/“+2(x) < C'l|f||H'>(x)- (30) 
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Such estimates also hold in Sobolev norms for every 1 < p < oo. Using a lift operator, if 
the boundary is smooth and problem (29) is replaced by a boundary value problem with f = 0 
and u|gx = g G (dX), one may obtain a similar inequality as (30) upon replacing ||f 






by ||g| 

Additionally, we will use below that when the elasticity tensor C is smooth and uniformly 
pointwise stable, interior regularity arguments as in [17, Theorem 1 p309] (in the case of scalar 
elliptic PDEs) translate into the following: if u solves (29), then for V CC U C X, there exists 
a constant C'i([/, V, C) such that 


e\\m{v) < C'i(I|u||i, 2([7) + ||f||i2(f/)) 



1 


(^diUj -\- djUi') 


(31) 


Finally, we note that for any C A, the functionals Xi : [C^(n, S' 3 (M))]® —>■ C®(n) and 
X 2 : [C\n, S 3 {R))f+^ C0(O), dehned by 

{ 6 (^)}i <,<6 ^ Ai(e(b,... ,e( 6 )) = det(e(b,... ,^( 6 )), ( 32 ) 

{eb)}i<,.<6+^ ^ M{M') : (33) 

M'CM,#M'=20 


are continuous as polynomials of the entries of the strain tensors and their hrst partial derivatives. 


5.2.2 Near-constant elasticity tensors 

We now use the constructions from the constant coefficient problem in Sec. 5.1 as well as the 
continuity of displacement fields with respect to the elasticity coefficients in appropriate norms, 
to establish reconstructibility of near-constant tensors. 

Proof of Theorem 2.7(i). For two displacement fields u, u' solutions of (1) with respective elas¬ 
ticity tensors C, C (both uniformly pointwise stable and at least C^-smooth) and the same 
boundary condition g G 772 (9A), the difference u — u' satisfies the following PDF 

V • (CV(u - u')) = V • ((C" - C)u') (A), (u-u')| 9 x = 0, 

so that, using (30) with s = 2, we obtain that 

Ik “ ^'Wmix) < K\\C - C"||c3(x)||u'||//3(x) < K'\\C - C"||c3(j)f)||g||^5^^^^, 

for some constants K,K'. By Sobolev inequality H^{X) —)■ C^’ 2 (A) C C^(A), we arrive at the 
estimate 


e-c'||ci(^)<A||C-C'||c3(x). 


(34) 
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Now pick C" = Co a tensor with constant coefficients and construct solutions sat¬ 
isfying hypotheses 2.2, i.e. ..., and \ ■ ■ ■, are bounded away from zero 

throughout X. We then see that combining the continuity of Fi,X 2 in (32)-(33) with estimate 
(34) applied to ..., shows that, if an elasticity tensor C is close enough to Cq in 

C^(X)-norm, then upon defining solutions of V • (CVu^-^^) = 0 with boundary condi¬ 
tion the difference — ec/^||ci(x) should be so small that 

and F 2 {€.^^\ • • •, remain bounded away from zero throughout X, so that C can be ex- 

plicitely reconstructed from these strain tensors. Theorem 2.7(i) is proved. □ 

5.2.3 Runge approximation 

We say that a differential operator C satisfies the Runge approximation property on X if for 
every compact fl C X, every solution v of Cv = 0 (fl) can be approximated in L^(fl) by solutions 
u of Cu = 0 (X). Other approximation topologies can be used (e.g. in [31, Theorem 

3.3]), in fact here we will use a approximation. The Runge approximation property is closely 
related to the unique continuation property, see e.g. [25]. 

Unlike for scalar elliptic PDE’s, the literature on unique continuation results for elliptic 
systems is sparse. In the context of elasticity, it has been proved for the Lame system (and 
more generally for elliptic systems with an iterated Laplacian as diagonal part) using doubling 
inequalities in [1], and for the transversely isotropic case in [31]. 

When the Runge approximation is available, following ideas from the first and third author 
presented in [13] and generalized in other contexts in e.g. [29, 9], one becomes able to fulfill 
the hypotheses of reconstructibility 2.2 by constructing local solutions satisfying appropriate 
qualitative behavior. This is what we now do in the proof of Theorem 2.7(ii). 

Proof of Theorem 2.7(ii). We decompose the proof in three steps. 

Step 1. Local solutions ^vith constant coefficients: Fix xq G X and Rsr = B^ri^o) C X 
a ball of radius 3r (r tuned hereafter) centered at xq, and denote Cq := C'(xo). Following 
the approach in section 5.1, we first construct solutions to the problem with frozen constant 
coefficients, i.e. such solutions, which we denote (each with strain tensor solve 

V • (CoeJ/'^) = 0 (M^), l<i <6 + W (35) 

This family is constructed so that Hypotheses 2.2.A-B are satisfied throughout B^r (in fact, they 
fulhll these hypotheses globally). 

Step 2. Local solutions with varying coefficients: From solutions we con¬ 
struct a second family of solutions (each with strain tensor via the following 
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equation 


V-(Ce(^)) = 0 {Bsr), l<j<6 + N. (36) 

The difference of both solutions satishes, for 1 < j < 6 + iV, 

V-(C'(e(^)-e(-^))) = V.((C7o-C)6(^'^) (B^r), = 0, (37) 

where the right-hand side is smooth since C is assumed to be smooth and so is Then 

using estimate (30) in high enough Sobolev norms combined with Sobolev inequalities [17, Th. 
6 p270], we arrive at estimates of the form 


u, 






< Cjlu 


( 1 ) _ 

T 


U, 


( 1 ) 


||H4(B3.)<C''||V-((Co-C)e(^'^)||H2(B3.) 


so that 


lim max -i = 0 . 

'■-s-o l<j<6+N " ' 


(38) 


Step 3. Runge approximation (control from the boundary dX): Assume r has been 
hxed at this stage. By virtue of the Runge approximation property, for every <5 > 0 and 
1 < j < 6 -b A^, there exists G H^(dX) such that 

11 ^ 2 ( 53 ^) < 5, where solves (1) with jax = (39) 

Since C is smooth and V • {C{e^p — = 0 thoughout B^r, estimate (31) with f = 0 implies 

that 

Il4^^ -4^^11^1(52.) ^ -uW||i 2 (B 3 ^). 

Moreover, using that V■{Cdpq{e^^^'^— er'^)) = —'V-{{dpqC){e\^'^ — ei^'^)) for every I < p < q < 3, and 
assuming that C is at least of class C^, we can hnd a constant 6*2 (r, C) (we drop the dependency 
on C below) such that 

Il4^^ - ^r^\\H\Br) < <^2(r',C)||e^^^ - e^^\\Hl[B2r) ^ <^3(0l|ui^^ - uW||i2(B3^) < C3(r)5. 

By Sobolev inequality —?■ C ^’2 C C^, we deduce that 

ll4'^-4'4lci(5.)<C4(r)5, l<j<6 + N. 


Since r is hxed at this stage, we deduce that 


lim max -i = 0 . 

5^0 i<j<6+N" ^ ^ 


( 40 ) 
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Completion of the argument: Theorem 2.5 states that C is reconstructible if both func¬ 
tionals J-i^J -2 are bounded away from zero for some properly chosen solutions. Fixing again 
xq € X and Br, as above, Step 1 established that ..., and • • •, were 

bounded away from zero over Br, where the Cq'^^’s were strain tensors associated with a problem 
with frozen coefficients. Due to limits (38) and (40), there exists a small r > 0, then a small 5 > 0 
and local solutions (u^^^ , • • •, such that maxi<j< 6 +Ar llci(Sr(a;o)) small that, 

by continuity of and B 2 mentioned above, 7u(e^^\ ..., e^^^) and B 2 {€^g\ ■ ■ ■, remain 

uniformly bounded away from zero over Br- Hypotheses 2.2 are thus satisfied over Br by the 
family which is controlled by boundary conditions. C is thus reconstructible over Br 

with a local stability as in Theorem 2.5. The proof of Theorem 2.7(ii) is complete. □ 
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